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We establish a general relation between the diagonal correlator of eigenvectors and the spectral 
. . . Green's function for non-hermitian random-matrix models in the large- Ai" limit. We apply this 

O^ ' result to a number of non-hermitian random-matrix models and show that the outcome is in good 

0^ ' agreement with numerical results. 
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I. INTRODUCTION 



Recently, a number of new results for non-hermitian random-matrix (NHRM) ensembles were obtained (see e.g. 
|0-p[ and references therein), reflecting on the rapidly growing interest in properties of NHRM in several areas of 
physics. In a recent paper |^, Chalker and Mehlig have pointed out the existence of remarkable correlations between 
left and right eigenvectors associated with pairs of eigenvalues lying close in the complex plane. Such effects do not 
*y-s ' exist for hermitian (more generally normal) random-matrix models, since in this case the left and right eigenvectors 
f\l , coincide. Some observables related to the eigenvector properties in non-hermitian random-matrix models have been 
f^ ' introduced and studied numerically in ^. Chalker and Mehlig Q obtained analytical formulae (in the large- A^ 
0^ ' limit, where N is the size of the NHRM) for correlations between left and right eigenvectors in the case of Ginibre's 
^^ ■ ensemble [pi. However an efficient calculational scheme for the simplest (and perhaps physically more transparent) 
one-point function 0(z) was lacking. In this paper we prove a simple relation stating that the correlator between left 
and right eigenvectors corresponding to the same eigenvalue is exactly equal, in the large- iV limit, to the square of the 
off-diagonal one-point Green's function |4| for non-hermitian eigenvalues. The latter is readily calculable for a wide 
n^ : variety of NHRM ensembles. We illustrate our observation in a number of NHRM models and show that it agrees 
(— I with numerical calculations. 

Q The present results can be used to study the interplay between reorganization of the left and right eigenvectors 

O ' and structural changes in the complex spectrum. In the simplest non-hermitian model - Ginibre's ensemble (or 

^ [ generalizations thereof 0,11^1) - the density of complex eigenvalues is constant, filling uniformly the circle (ellipse) 

in the complex plane. This model and its variants do not have external parameters, which could induce structural 

changes ("phase changes") in the average eigenvalue distribution e.g. the changes from simply to multiple-connected 

^ , domain of eigenvalues. In order to study these more complex phenomenae one has to consider e.g. the model for open 

■ - ' chaotic scattering, where at large couplings (strong dissipation) the eigenvalue distribution splits into two disconnected 

islands, reflecting on the separation of time scales |[ll],12l. One island corresponds to short-lived resonances, while 

the other to long-lived (almost classical) trapped states |13|]. It is known in this case, that the reorganization of the 

eigenvalues is followed by some reorganization of the eigenvectors m. In particular, the norm of the right states is 

sensitive to the distribution of resonances . 

In section 2, we introduce the notations and summarize the main results for the eigenvector correlators in the case 
of Ginibre's ensemble g|, as established recently by Chalker and Mehlig Q. In section 3, we present our main result 
and apply it to several NHRM models with direct comparison to numerical results. Our analysis relies on novel 
techniques for NHRM models discussed by some of us Q. A summary of our results is given in section 4, and a 
number of technical details can be found in the Appendices. 

II. GINIBRE'S MATRIX MODEL 

Ginibre [g| has introduced a Gaussian ensemble of general complex matrices N x N, i.e. matrices distributed with 
the probability 



P{M)dM - exp{-NTTMM'')dM 



(1) 



giving non-vanishing cumulants {MabMat) = 1/-/V. The eigenvalues are uniformly distributed on a unit disc centered 
at the origin of the complex plane. In Appendix B, we provide a short derivation of this result and others using 
(matrix- valued) Blue's functions |H,p[. Since the matrices are complex (non-hcrmitian), there exists a bi-orthogonal 
set of right (i?) and left (L) eigenvectors, so that 



M=Y,>^a\Ra><La\ , 

a 

Alt^^Afc|L6><i?b| , 

b 

where < La\Rb >= Sab- Chalker and Mehlig |q| have studied the following eigenvector correlators: 

oiz) = (^j2o--^('-^-n ' 



(2) 



0{z,w) 



a^b I 



(3) 



where Oab =< La\L}j >< Rb\Ra >■ The main results of their paper are exact expressions for 0{z,w) and for 0{z) in 
the case of Ginibre's ensemble. For N ^ l,\z — w] ^ 0, where z, w are lying within the unit circle 

1 1 — zw 



0{z,w) ~ — - 



For \z 



0, 



0{z,w)r 



N^1^\Z\ 



.[i-(i+ic.nexp(-ic.n] 



where Z = {z + w)/2 and u = \/N{z — w). The diagonal correlator reads 

0(z) = -(l-|zp). 



(4) 



(5) 



(6) 



Whereas the calculation of the "wide" eigenvector correlators 0{z, w) (like (|J)) in many cases of NHRM is straightfor- 
ward, using e.g. non-hcrmitian diagrammatics H , the calculation of the "close" eigenvector correlators 0{z — w)micro 
and the diagonal correlators 0{z) is technically involved. In practice, it seems that explicit calculations in the micro- 
scopic limit are possible only in the cases when the spectrum of NHRM possesses an azimuthal (rotational) symmetry. 
In the next section we provide a general formula for the diagonal eigenvector correlator 0{z) in terms of the spectral 
one-point Green's function. 
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FIG. 1. Eigenvector correlator 0{z = 0-1- iy), generated from an elliptic ensemble of 5-10^ matrices of size 20 (crosses), an 
ensemble of 10^ matrices of size 40 (stars) and size 100 (boxes) versus the analytical correlator —QqqGqq, for r — 0.5 (left) and 
T = 0.7 (right). 



III. NEW FORMULA AND RESULTS 

The main result of this paper reads 

N 

oiz) = — g,,-g,,. (7) 

TT 

Here Qqq and Qqq are ofF-diagonal elements of the generalized (2 by 2) spectral Green's function Q |Q] 

g = ( f,'^ f,'^) . (8) 

The elements Qab are defined as traces Qab = j^'^'^NQab of the N x N blocks of the generalized resolvent Q, 

Gqq Gqq\ / f Z- M IeIat 



where Ijv is the TV-dimensional identity matrix and < ... > denotes the averaging over the pertinent ensemble of 
random matrices Ai. In comparison to the original work B, we have chosen here purely imaginary infinitesimal 
values in the off-diagonal block. This way guarantees that the mathematical operations performed in the proof of (0) 
(cf. Appendix A) are well-defined. 

The spectral density follows from Gauss law [Q, 

V{z,z) = -dg gqq{z,z) (10) 

TT 

which is the distribution of eigenvalues of Ai. For hermitian Ai, ( |lO|) can be nonzero only on the real axis. As 
e ^ 0, the block-structure decouples, and we are left with the original resolvent. For z -^ +iO, the latter is just a 
measurement of the real eigenvalue distribution. 

For non-hermitian Ai, as e — > 0, the block structure does not decouple, leading to a non-holomorphic resolvent for 
certain two-dimensional domains on the z-plane. For more technical details we refer to the original work H, or recent 
reviews 0,^ . Similar constructions have been proposed recently in |^,^ . 

The r.h.s. of the relation (^ is usually given by a simple analytical formula. Technically, the most efficient 
way of calculating the off-diagonal components of the Green's functions is to use the generalized Blue's function 
technique [H,p|. In Appendix B we provide a pedagogical derivation of some of the results below, for others we refer 
to the original papers. 

For Ginibre's ensemble we immediately get (cf. Appendix B) 

Oiz)G^„ibre = \V^^^Tf = —(1 - \zf) (11) 

TT TT 

in agreement with Chalker and Mehlig g|. For the elliptic ensemble ^,|l^, using the off-diagonal elements from 
Appendix B, the diagonal correlator reads 

0{z)EmpHc = ^ .^^^2)2 [(1 - ^')' - (1 + ^')l-^l' + 2^1^602] . (12) 

In Fig. n, we present numerical results generated from computer simulation of eigenvectors for the elliptic ensemble 
with T = 0.5 and r = 0.7 and different size of the matrix, N, versus the analytical prediction (|l2|). The results are 
satisfactory. The figure provides also a rough estimation of finite-size effects. We note that numerical simulations of 
eigenvector correlations are time consuming, hence the utility of the result (|^), where the r.h.s is straightforward to 
calculate. 

The models considered above have constant density of complex eigenvalues, and the domain of eigenvalues is simply 
connected (circle or ellipse). Below, we consider a toy-model, where the domain could split into two disconnected 
domains at some critical value of the external parameter. Also the distribution of eigenvalues is non-uniform. The 
simplest non-hermitian model of this kind is Ginibre's random Hamiltonian plus a two level deterministic Hamiltonian, 
with A^/2 levels a and N/2 levels — a.R Using the addition law for the generalized Blue's function, we easily obtain all 
the components of the matrix valued Green's function Q (cf. Appendix B). 



^To our knowledge, this model was first considered by Feinberg and Zee, using their hermitization method H. 



Using the relation (^ we predict 

0{z)g+d 



N 

■K 



-(l + v/l + 4a2(z + z)2) 



(13) 
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FIG. 2. "Ginibre+detenninistic" eigenvector correlator 0{z = + iy), generated from an ensemble of lO"" matrices of size 60 
(crosses) versus the analytical correlator —GqqGqq, at a = 0.5 (left), a = 1 (center), and the correlator 0{z = a; + iO) at a = 2 
(right). 

We note that the value a = 1 is the critical one, when the single island of eigenvalues splits into two. Numerical 
simulations for this ensemble are shown in Fig. ||. The solid line is the analytical result, the crosses are the numerical 
results calculated for the line z = + iy. The diagonal eigenvector correlators 0{z) inside the islands follow the 
distribution determined by the off-diagonal components for the spectral Green's function. We also note that the 
eigenvector correlator follows precisely the boundary-shape of the eigenvalues. This is expected, since the condition 
\Qqq\ = determines M the regions in the complex plane separating the holomorphic and non-holomorphic components 
of the spectral Green's function q 

Finally, we consider the case of open chaotic scattering. As an illustration we choose the classical result of Haake, 
Sommers and coworkers ||l2|, based on Mahaux-Weidenmiiller |ll) microscopic picture for nuclear reactions. In brief, 
the model is generically described by a non-hermitian Hamiltonian of the form 



H - igVV^ , 



(14) 



where _ff is a random Gaussian (orthogonal) N x N matrix, while V is an N x M random matrix |T2|. Here N 
is a number of discrete states and M < A'^ is a number of the continua. The model was solved [|l2|in the limit 
N — + oo, M -^ oo, m = AI/N fixed. Using the results from Blue's function techniques (cf. Appendix B) we predict 
the analytical behavior for the correlator 0{z) for this model as 



0(z) 



N 



1 



1 - gy 



4 VI -52/ 



ra 

y 



(15) 



with z — X -\- iy. In Fig. we compare this result to a numerically generated ensemble of eigenvectors. We would like 
to note large finite size effects at the edges of the islands. 

Since in this model a change in the parameters g and m can generate structural changes in the distribution of 
eigenvalues, it is interesting to study whether the splitting of the "islands" is accompanied by some distinct behavior 
of the eigenvector correlators. 



■^The shape of the eigenvalue domains for NHRM models can be inferred from associated hermitian models using conformal 
mapping |4| . This points to yet another relationship between the eigenvalues and eigenvectors of hermitian and non-hermitian 
models. 
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FIG. 3. Scattering model eigenvector correlator 0{z = + iy) generated from an ensemble of 10^ matrices of size 60 (crosses) 
and of size 120 (boxes) with m = 0.25, versus the analytical correlator —QqqQqq, at different values of the coupling constant g. 
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FIG. 4. Three dimensional plots of the analytical results for the eigenvector correlator 0{z), normalized by the spectral 
density v[z, z), versus x and y for different couplings g un the case of open chaotic scattering. 

In Fig. ^ we plot the analytical result (|5|), normalized by the spectral density i'{z,z) (cf. Appendix B), as a 
function of x and y, where z — x + iy. For the case considered here (m = 0.25) the value of the critical coupling is 



g = 4.44.0 We observe a strong reorganization in the distribution of the average "norm" of eigenvectors in the vicinity 
of the critical couphng. 

We have not taken up in this paper the issue of wide eigenvector correlators 0{z,w). For the cases considered 
here, these correlators are readily constructed using the NHRM diagrammatic Bethe-Salpeter equations Q, as noted 
by Chalker and Mehlig H. In most cases, however, the resulting final formulae are rather lengthy. We note that 
the knowledge of 0{z,w) in the regime z — w ^ 0(1) does not suffice to determine 0{z) through the sum rules |g] 
originating from the bi-orthogonality of the left/right eigenvectors. Hence the relevance of the present investigation. 

IV. SUMMARY 

We have established a novel relation between the diagonal correlator of eigenvectors and the off-diagonal elements 
of one-point spectral Green's function for general ensembles of NHRM models. We have appplied this result to a 
number of NHRM models and checked its validity against numerically generated results. Our observation accounts 
for part of the eigenvector correlations established recently by a number of authors |g,0| . Our result generalizes to 
non-hermitian ensembles with real, complex or quaternionic components, as well as non-hermitian ensembles with 
additional symmetry (e.g. chiral NHRM). In this last case, however, non-hermitian diagrammatic techniques have to 
be used instead of the versatile method of Blue's functions. 

Finally, we point at the possibility of relating the eigenvector correlators to the eigenvalue correlators in the 
microscopic limit. This issue and others will be discussed elsewhere. 
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APPENDIX A: PROOF OF EQ. (7) 

We prove first the following representation for 0{z) 



lim ( TrT ---, — rr ^ • Tr^ — r-- — ^ ) = TriV 0{z) . (Al) 

e-*o\ {z~M){z-M^) + e^ {z-M^){z-M)+e^ / ^' ^ ^ 

Below we use Ti A — logdetA, and the fact that M can be diagonalized by a non-unitary transformation U, 
U~^MU = diag (Ai, . . . , A^v). Introducing V = IJ-^U-'^^ we get 

Using this representation, we note that 0(z) is zero in the limiting procedure, unless z is close to any (say Ai) of the 
eigenvalues of A^. If this happens, we use the parameterization, z — Xi = eu exp{i(j)) with u ^ 0{l). Of course we get 
additional (similar) contributions when z is close to the other eigen valu es A.;. For notational simplicity we will now 
consider only the case of Ai, adding the remaining contributions in ([A7). 



Using Laplace's expansion for the first row of the determinant in (A2) we arrive in the e — > limit 



-9,logdet[---l = i-TuT- , (A3) 

^ e i-hu ViidctWi 

with W\ being the matrix (z — Ai)Vifc(z — Afe) without the first row and column and (M^2)i+i fc+i = (W^i)ifc, (1^2)1 k = 
(z - Afe)Fu, [W2)^^ = {z- Xk)V,i and (1^2)11 ^Vn. 



^Generally, gl,, = (1 - ^)-« |I|. 



From Laplace's expansion 

detW'i^dety' Yl {z - X^){z ~ X,) , 

i=2,N 

det W2 = det y \[ {z - X,){z ~ Xi) (A4) 

i=2.N 

the ratio of the determinants det Wi/ det W2 is simply (y~^)ii. In the above, V' is the minor of V left after crossing 
out the first row and first column. In this way we obtain: 

^9.1ogdet[...] = l ' , , (A5) 



The determinant corresponding to the second trace in ( Al), with (z — A4) and {z — A4^) interchanged is obtained by 
substituting V <-> V~^ and hence is given by exactly the same formula. 
Therefore 



(^ttS^^^z - Xi)Vii{V-')ii) . (A6) 



where we used the representation for complex Dirac delta nS^^^z) — lime^o e^/(e^ + kP)^- 
In this way we obtain the important formula 



i(9,logdet[...]9,logdet[...])-7r(^l^,,(y-i),,,52(z-A,)), (A7) 



where we have reinstated the sum over all eigenvalues. It remains to show that Vii{V ^)ii =< Ri\Ri >< Li\Li >. 
Using any orthonormal basis {ci} and the decomposition (0), we see that the linear transformation U satisfying 

U~^MU = Y,>^k\ek><ek\ (A8) 

k 

can be written explicitly as 

U = Y,\Rk><ek\ , 

k 

y-i = [/t[/^^|efc><i?fc|i?„><e„|. (A9) 

k,n 

From the last equation we infer 

lA-i =< e,|y-i|e, >=< i?,|i?, > . (AlO) 

Similar reasoning leads to Vu =< Li\Li >. 
This yields 

Va, logdet[---] a, log det [•••]) =7r/^ < L,\L, >< R,\R, > S^^\z ^ X,)\ = ttNO{z) (All) 



which completes the proof of (Al) 



We would like to stress that till now we have not used the large- A'' arguments in deriving this formula. Therefore 



(Al) is more general than the relation (|7|), and could be used as a starting point for a systematic study of finite-size 
effects or microscopic limit. 

To complete the proof of relation (|^), we observe that, in the large- A^ limit, we could use the factorization theorem, 
such that the l.h.s. of ( |AlJ}^ splits into the product of averages < Tr[...] >< Tr[...] >. 

From the definitions (H,0), the off-diagonal Green's functions have the following expression: 



Gqq = ( -J lim — Tr 



ON {z - M){z ~ M^f) + e^ 
5 7V \z - Ml){z ^ M) + e^ 



S,, = ( -* lim -Tr ,^ ^^,,,^ ^^, , ^, ) (A12) 



Hence 



TT 

TV 



-0(z) = -g^gGgq (A13) 



which completes the proof of equation (^ . 



APPENDIX B: GENERALIZED BLUE'S FUNCTION 

The generahzed Blue's function ^^ is a 2 x 2 matrix valued function defined by 

B[g{Z)]=Z=(^^^ I') , (Bl) 

where G was defined in Section 2 and e will be eventually set to zero. This is equivalent to the definition in terms of 
the self-energy matrix 

Big) = i: + g-\ (B2) 

where S is a 2 x 2 self-energy matrix expressed as a function of a matrix-valued Green's function. The addition law 
for generalized Blue's functions reads 

Z ^ B.iG) + 82(0) - g~' . (B3) 

in analogy to the original construction by Zee [|l7| for hermitian matrices. 

1. Ginibre's ensemble 

Ginibre's ensemble [g| could be viewed as a sum of hermitian and anti-hermitian Gaussian ensembles, with the 
original width suppressed by \/2 in relation to the original width of the complex gaussian ensemble. The generalized 
Blue's function for the hermitian part is simply ^ 

Br{A) = ^A + A-\ (B4) 

The generalized Blue's function for anti-hermitian part is B 

B,r{A)^^A + A-\ (B5) 

where we used the notation 

The factor 1/2 comes from the normalization of the width, and the extra signs reflect on the anti-hermitian correlations 
of the matrix elements. The addition law reads now (in the matrix form) 

z = g-' + ^[g + g]. (bt) 

The nontrivial (non-holomorphic solution) reads 



»^l/FF^^^t^l^ (B8) 



The domain of eigenvalues is determined by the condition Qqq = 0, for which the block structure decouples leading 
to holomorphic and anti-holomorphic copies. For this ensemble, this is simply a circle |zp = 1. Inside the circle, 
Qqq{z, z) = z (upper left corner of G). The constant density of eigenvalues follows from Gauss law (|o|). Outside the 
circle, the second (holomorphic) solution of (B7) is valid, giving G{z) = 1/z. This reproduces the salient features of 
Ginibre's ensemble. 

A straightforward generalization jj] using the measure ||9||l0|] 



P{M)dM - exp 



N 



l-T- 



:Tr(AlAlt _ tRcMM) dM 



(B9) 



leads to the results for the elliptic ensemble with 

Z — TZ 



Gqq = 



1 



Gqq = G, 



GqqGq 



qq , 



2tRcz^ - (1 - r2)2 



(1-r^ 



(BIO) 



Note that the measure (89) leads to non- vanishing cumulants (MabMab) — '^/N and (MabMba) — t /N . In particular, 
T = —1 corresponds to anti-hermitian matrices, explaining the flips of the signs in the tilted variables above. 



2. Tavo level deterministic Hamiltonian plus Ginibre's ensemble 



Since the deterministic hermitian Green's function for the two-level Hamiltonian is 

Gd{z) ^ ^ ^ ^ 



2 \z — a z + a^ 
the corresponding generalized Green's function for this Hamiltonian reads 

GoiZ) = \[{Z- al2)"' + {Z + aU)-^] 



with I2 denoting the two-dimensional identity matrix. Substituting in (B12) Z —^ Bd{G) we obtain 



e = 2 [i^oiG) - aU)-^ + {Bd{G) + al2)~'] 
The addition law for generalized Blue's functions reads 

B{A) = Bd{A) + Bg{A) - A-^ 



(Bll) 



(B12) 



(B13) 



(B14) 



where the generalized Blue's function for Ginibre's ensemble was constructed above, i.e. Bg{A) = .A ^ + 1/2{A + A). 
Substituting in (B14) A -^ G{Z) we infer the relation 



Z^Bd{G) + ^{G + G) . 



From (B13) and (B15) we get the final equation 



1 



{Z-^{G + G)- al2)-^ + {Z- ^{G + G) + al2)-^ 



(B15) 



(B16) 



Solving this matrix equation we arrive at 

GqqGqq - 



|zp +0"- -(1 + Vl + 4a2(z + z)2) , 



2a^{z + z) 



1 + a/1 + 4a2(z + z)2 ' 



Gqq ~ G, 



qq — Vqq ■ 



(B17) 



3. Open chaotic scattering 

Since the addition law for open chaotic scattering was formulated by us in previous publications 0] , we will be brief 
and refer to the original papers for details. The addition law reads 



z = m{i ^ rgy'^r + g + g-^ 



(B18) 



where m — M/N and T — diag (—«<?, 19). Solution of the matrix equation (B18) leads to 



Sqq 


X i 

= 2 + 2 


\ m q 

1 1 


.9 ' V ' 1 - 9V. 


Sqq 


x^ 1 
= T + 4 


( q ml 

1 


V1-.92/ ' y '9 


Qqq 


^Sqq 





where z = x + iy. Gauss law leads to the spectral density | [l^ 



A'Kv{x^ y) = div E = 1 -\ — ^ ^ 



1 



1 - 9y 



y^ (i~5y)^ 



(B19) 



(B20) 



where E^ = 2Re gqq and Ey = — 2Im ggg. For completeness we mention that condition ggg = reproduces the results 
by Haake, Sommers and coworkers |12] for the boundary of eigenvalues in open chaotic scattering. 
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